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Equation of state of baryon rich quark matter is studied within the 
SU(3) Nambu-Jona-Lasinio model with flavour mixing interaction. Possi- 
ble bound states (strangelets) and chiral phase transitions in this matter are 
investigated at various values of strangeness fraction r s . The model predic- 
tions are very sensitive to the ratio of vector and scalar coupling constants, 
£ = Gy/Gs- At £ = 0.5 and zero temperature the maximum binding energy 
(about 15 MeV per baryon) takes place at r s ~ 0.4. Such strangelets are 
negatively charged and have typical life times ~ 10 s. The calculations 
are carried out also at finite temperatures. They show that bound states 
exist up to temperatures of about 15 MeV. The model predicts a first order 
chiral phase transition at finite baryon densities. The parameters of this 
phase transition are calculated as a function of r s . 



I. INTRODUCTION 

Almost 30 year ago A.B. Migdal has put forward a brilliant idea of pion condensation 
in nuclear matter [I]]. Soon it was realized that this phenomenon may lead to the existence 
of density isomers ||. At about this time Lee and Wick proposed another mechanism 
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leading to the appearance of an abnormal nuclear state || . It is related to the restoration 
of chiral symmetry at high baryon density. These ideas have motivated nuclear community 
to initiate new experimental programs aimed at producing hot and dense nuclear matter 
in energetic collisions of heavy nuclei. These experiments have started in Dubna and 
Berkeley and then continued in Brookhaven and CERN. Nowadays exciting expectations 
are associated with new ultrarelativistic heavy ion colliders, RHIC and LHC. 

The general goal of present and future experiments with ultrarelativistic heavy ions 
is to study the equation of state and dynamical properties of strongly interacting matter. 
Now the main interest lies in investigating the chiral and deconfinement phase transitions, 
predicted by QCD. The ultimate goal is to produce and study in the laboratory a new 
state of matter, the Quark-Gluon Plasma (QGP). This state of matter can be reached 
only at high temperatures or particle densities when elementary constituents, quarks and 
gluons are liberated from hadrons. 

Since the direct application of QCD at moderate temperatures and nonzero chemical 
potentials is not possible at present, more simple effective models respecting some basic 
symmetry properties of QCD are commonly used. One of the most popular models of 
this kind, which is dealing with constituent quarks and respects chiral symmetry, is the 
Nambu-Jona-Lasinio (NJL) model @||. In recent years this model has been widely used 
for describing hadron properties (see reviews 0,01), phase transitions in dense matter 
H| [13| and multiparticle bound states [Ll-|l? . 



In the previous papers [0,0 we have used the NJL model to study properties of 
the quark-ant iquark plasma out of chemical equilibrium. In fact, we considered a system 
with independent densities of quarks and antiquarks. We have found not only first order 
transitions but also deep bound states even in the baryon-free matter with equal densities 
of quarks and antiquarks. In the present paper the emphasis is put on investigating the 
possibility of bound states and phase transitions in equilibrated matter at various flavour 
compositions. In particular, we consider the possibility of bound states in quark matter 
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with a significant admixture of strange quarks, i.e. strangelets. Thermal properties of 
strange quark matter are also studied. 

The paper is organized as follows: in Sect. II a generalized NJL model including 
flavour-mixing terms is formulated in the mean-field approximation. The model predic- 
tions for strange matter and characteristics of its bound states at zero temperature are 
discussed in Sect. III. Effects of finite temperatures are considered in Sect. IV. Possible 
decay modes of new bound states are discussed in Sect. V. Main results of the present 
paper are summarized in Sect. VI. 



II. FORMULATION OF THE MODEL 



Below we use the SU(3)-flavour version of the NJL model suggested in Ref. [20]. The 
corresponding Lagrangian is written as (% — c = 1) 



j=0 

8 

-Gvj: 

j=0 



K [det/ (> (1 - 75) V) + det/ (> (1 + 7s ) $ 



Here ip is the column vector consisting of three single-flavour spinors ipf, f = u,d,s, 
Ai,...,Ag are the SU(3) Gell-Mann matrices in flavour space, Aq = 



'2/3 J, and 

m = diag(m 0tt , m 0( i, ttiqs) is the matrix of bare (current) quark masses. At m = this 
Lagrangian is invariant with respect to SUl(3) ® SUr(3) chiral transformations. The sec- 
ond and third terms in Eq. (|l|) correspond, respectively, to the scalar-pseudoscalar and 
vector-axial-vector 4-fermion interactions. The last 6-fermion interaction term breaks 
the Ua(1) symmetry and gives rise to the flavour mixing effects. 

In the mean-field approximation the Lagrangian ([!]) is reduced to 



G, 



3 



+ ^T,Pvf 2 + ^l[Psf, (2) 
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where Ip = <j) + i "foGvPvf an d 



Psf = <ip f ipf>, (3) 
pvf = < ipfjoipf > (4) 

are scalar and vector densities of quarks with flavour / . Angular brackets correspond to 
the quantum-statistical averaging. The constituent quark masses, to/, are determined by 
the coupled set of gap equations 

m f = wo/ - Gspsf + 2K Y[ Psf ■ (5) 

f'H 

The NJL model is an effective, non-renormalizable model. To regularize the divergent 
contribution of negative energy states of the Dirac sea, one must introduce an ultraviolet 
cut-off. Following common practice, we use the 3-momentum cut-off 0(A— p) in divergent 
integrals Q. The model parameters Too/, Gs, K, A can be fixed by reproducing the observed 
masses of ir, K , and rj' mesons as well as the pion decay constant f n . As shown in Ref. [BtJ, 
a reasonable fit is achieved with the following values: 

m 0u = m od = 5.5 MeV, m 0s = 140.7 MeV, (6) 
G s = 20.23 GeV" 2 , K = 155.9 GeV~ 5 , A = 0.6023 GeV. (7) 



Motivated by the discussions in in Refs. fl6|j2Tf , we choose the following value of the vector 
coupling constant^ 



G v = 0.5 G s = 10.12 GeV~ 2 



1 Here 6 (x) = ^ (1 + sgnx). 

2 See discussion of this question in Ref. [19] 
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Let us consider homogeneous, thermally (but not, in general, chemically) equilibrated 
quark-ant iquark matter at temperature T . Let ap t \ ( &j>,a) and ap~ x ( bp X ) be the de- 
struction and creation operators of a quark (an antiquark) in the state p, A , where p is 
the 3-momentum and A is the discrete quantum number denoting spin and flavour (color 



indices are suppressed). It can be shown [18] that quark and antiquark phase-space 
occupation numbers coincide with the Fermi-Dirac distribution functions: 



< cip X ap t \ >= ripf 
< b p,x h pXf >= npf 



( Epf ~ PRf \ , i 
exp — - ) + 1 



exp 



T 

E Pf ~ ~Pnf 
T 



-i 



-\ -i 



(9) 
(10) 



where Epj 



m 



f 



p 2 and [iRf, [i R f denote the reduced chemical potentials of quarks 



and antiquarks: 



P-Rf — — G v pvf , 
~pRf = 7 1 / + GvPvf ■ 
The explicit expression for the vector density can be written as 



Pv f = p f -p f , 



where 



p f = v 



d 3 p 



n 



(27r)3 -P/' Pt-»J J^y n Pf 

are, respectively, the number densities of quarks and antiquarks of flavour / and v = 
2N C = 6 is the spin-color degeneracy factor. The net baryon density is obviously defined 

as 



v 



d 3 p 



■ Til 



(11) 

(12) 
(13) 
(14) 
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PB 



oT,PVf 



(15) 



The physical vacuum (pf = J>j = 0) corresponds to the limit npf = npf = . 

In general the chemical potentials pf and JIj are independent variables. The assump- 
tion of chemical equilibrium with respect to creation and annihilation of qq pairs leads to 
the conditions 



H f = -fif, f = u,d,s. (16) 

These conditions automatically follow from the relations valid for any qq system chemi- 
cally equilibrated with respect to strong interactions: 

Hi = Biii B + Si ii s + Qi V-Q (i = u, d, s, u, d, s) . (17) 

Here three independent chemical potentials, /ie, fis, Hq are fixed by the net baryon num- 
ber B , strangeness S and electric charge Q of the system. 

In heavy-ion collisions at high bombarding energies the partonic matter created is 
characterized by B ~ and S ~ 0, therefore Jlf ~ /// ~ 0. On the other hand, 
strangelets are finite droplets with nonzero B and S . Using Eq. ([H]) one can see that in 
this case hb, lis 7^ and, therefore, the inequality fi s > ii Uj d should hold. This conclusion 
shows that strong fluctuations are needed for strangelet formation in high energy nuclear 
collisions. 

If strangelet life times are long enough, the equilibrium with respect to weak processes: 
s^u + e~+V e , u + e~^s + v e , s + u <->• u + d, (18) 
may be also achieved. Assuming that /i e = = one arrives at the following conditions 

Hu = (ld = Hs- (19) 

As will be shown below, constituent masses m s of strange quarks in (mechanically) stable 
strangelets exceed chemical potentials of it, d quarks. At ii s > m s and T — » the (3 
equilibrium conditions ( |19D can be realized only after all s quarks decay. 

Using conditions fll6|) in the limit of zero temperatures one can see that the density of 
antiquarks vanishes in baryon-rich chemically equilibrated matter, p — > . In the mean- 
field approximation the reduced chemical potential of quarks with flavour / coincides with 
their Fermi energy: 

fiRf = \Jm}+p* Ff , (20) 
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where 



v 



is the corresponding Fermi momentum. 

Within the NJL model the energy density and pressure of matter as well as the quark 
condensates psf contain divergent terms originating from the negative energy levels of the 
Dirac sea. As noted above, these terms are regularized by introducing the 3-momentum 
cutoff 9 (A — \p\) . Then the scalar density is expressed as 



Psf = " 



r d 3 p mt r _ , 



(22) 



The energy density corresponding to the Lagrangian (0) can be written pi as 



e = e K + e D + e s + e v + e FM + e • (23) 



This expression includes: 
the "kinetic" term 



e K = v E / 7^4 E pf (n pf + np f ) , (24) 



/ 

the "Dirac sea" term 



(2k) 



e D = ~vY. J -0j- 3 E PfH^-p)i (25) 



the scalar interaction term 



*5 = ^ (26) 



the vector interaction term 



ev = ^T,Pvf (27) 



and the flavour mixing term 



e FM = -AK\{ Psf . (28) 
/ 



A constant e is introduced in Eqs. (^) in order to set the energy density of the physical 
vacuum equal to zero. This constant can be expressed through the vacuum values of 
constituent masses, m™ c , and quark condensates, pjgf. These values are obtained by 
selfconsistently solving the gap equations (]5|) in vacuum, i.e. at ripf = npf = . 

Explicit analytic formulae for the energy density and the gap equations may be ob- 



tained in the case of zero temperature. At T — > one has [18 



ex + e D 



v 



■ 7T Z 



E 



J>Ff) V A 



(29) 



E 



PFf 



PFf, 



A 



(30) 



where 



= A J dttWt 2 +X 2 = (l + VTT 



— x\ i + VT+x 2 

x z in 

2 x 



(31) 



For a system with independent chemical potentials for quarks (fif) and antiquarks 



(p,f) one can use the thermodynamic identity for the pressure of qq matter 

p = E^/P/ + V-fPf) -e + sT, 
f 



(32) 



where s is the entropy density 
d 3 p 



s = ~ v E 



(2vr) : 



n pf Innpf + (1 - n pf ) In (1 - n pf ) + n pf -> n p/ 



(33) 



As discussed in Ref. [Q, bound states of gg matter and first order phase transitions in 
qq matter are possible if its equation of state P = P(/j,f,~pf, T) contains regions with nega- 
tive pressure or isothermal compressibility, respectively. At T = the state of mechanical 
equilibrium with vacuum (P = 0) corresponds to minimum of the energy per particle, 
e = E / J2f(Nf + Nf) — ej J2f(Pf + Pf) > where is the number of quarks (antiquarks) 
of flavour / . In the case of pure quark matter Nf = 0, its baryon number B = A//3 
and its energy per baryon E/B = 3e . 
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To characterize the flavour composition we introduce the strangeness fraction param- 
eter 

Below only the isospin-symmetric mixtures where N u = Nd and iV s = iVj will be consid- 



ered. It can be shown [19] that in the dilute limit, when all single flavour densities 

are small, e tends to the sum of the constituent quark and antiquark masses in vacuum, 

weighted according to r s : 

m™ c {r s ) = (1 - r s ) m™ + r s (35) 

In the case of chemically equilibrated matter at T = and fixed r s one has 

e(p B ^0,r s ) = m™ c (r s ) (36) 

A bound multiparticle state exists if there is a nontrivial minimum of e as function of ps 
and the binding energy (BE) per baryon is positive: 



BE = 3 



m V cT(r s ) - e min (r s ; 



> . (37) 



III. STRANGE QUARK MATTER AT ZERO TEMPERATURE 

Let us consider first quark matter with nonzero net baryon density at T = 0. In 
the chemically equilibrated system the density of valence antiquarks will be zero for each 
flavour (jif = 0). Fig. 1 shows the energy per quark as a function of baryon density, 

Pb — \ J2Pf ■ Different curves correspond to different r s , which in this case is the relative 

/ 

concentration of strange quarks. In accordance with Eq. (|3^), at p# —>■ the energy per 
quark tends to the corresponding vacuum mass. With growing density both the attractive 
scalar and repulsive vector interactions contribute to e (see the discussion of this question 
in Ref. @). 
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It is interesting that at r s < 0.7 the attractive interaction is strong enough to produce 
a nontrivial local minimum at a finite p#. In the pure u, d matter (r s = 0) this minimum is 
unbound by about 20 MeV as compared to the vacuum masses of u and d quarks. On the 
other hand, it is located at a baryon density of about 1.8 p , which is surprisingly close to 
the saturation density of normal nuclear matter. Of course, the location of this minimum 
depends on the model parameters. Nevertheless, one can speculate that nucleon-like 3 
quark correlations, not considered in the mean-field approach, will turn this state into 
the correct nuclear ground state. 

When r s grows from to about 0.4, the local minimum is getting more pronounced and 
the corresponding baryon density increases to about 3.2 po. At larger r s , the minimum 
again becomes more shallow and disappears completely at r s ~ 0.7. At 0.2 < r s < 0.6 
the minima correspond to the true bound states, i.e. the energy per quark is lower than 
the respective vacuum mass. But in all cases these bound states are rather shallow: even 
the most strongly bound state at r s ~ 0.4 is bound only by about 15 MeV per baryon. 
Nevertheless, the appearance of local minima signifies the possibility for finite droplets 
to be in mechanical equilibrium with the vacuum at P = 0. It is natural to identify 
such droplets with strangelets, which are hypothetical objects made of light and strange 
quarks f2|-g7]]. 

It should be emphasized here that /3-equilibrium is not required in the present ap- 
proach (see the discussion below). That is why our most bound strangelets are pre- 
dicted to be more rich in strange quarks (r s > 1/3) than in the approaches assuming 
/^-equilibrium p2] , [2"3p5| , which give r s < 1/3. As a result, these strangelets will be neg- 



atively charged^. Indeed, the ratio of the charge Q to the baryon number B is expressed 
through r s as 

£ 2£. 1 A, 1 ft 1 (1 _ 3t , ). (38) 
B 3 p B 3 p B Sp B 2 



^Negatively-charged strangelets have been also considered in Refs. p4lE 
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For r s ~ 0.4 this gives Q/B ~ —0.1. In light of recent discussions (see e.g. Ref. [27] ) 
concerning possible dangerous scenarios of the negatively-charged strangelet production 
at RHIC, we should emphasize that the strangelets predicted here are not absolutely 
boundF], i.e. their energy per baryon is higher than that for the normal nuclear matter. 
Hence, the spontaneous conversion of normal nuclear matter to strange quark matter is 
energetically not possible. 

Fig. 2 shows the constituent masses of u and s quarks as functions of baryon density. 
The dropping masses manifest a clear tendency to the restoration of chiral symmetry at 
high densities. The dots indicate the masses at the local minima in the respective energies 
per baryon shown in Fig. 1. Note that the stronger is reduction of constituent masses the 
deeper are the corresponding bound states. For the metastable state at r s = 0, which is 
a candidate for the nuclear ground state, the masses of u and s quarks are equal to 0.3 
and 0.9 of their vacuum values respectively. For the most bound state at r s ~ 0.4 the 
respective mass ratios are reduced to 0.15 and 0.6. 

The behaviour of the u and s chemical potentials is shown in Fig. 3. At r s ^ 0, 1, 
due to the contribution of the vector interaction (see Eqs. (|TTD, (p0|) ), /i u and /i s grow 
practically linearly at large baryon densities. In accordance with a previous discussion, 
one can see that the conditions /i s > \i u ^ hold at baryon densities corresponding to bound 
states of strange matter. 

The properties of the multiparticle bound states are summarized in Figs. 4-5. Fig. 4 
shows the binding energy per baryon, Eq. fl3~T|). The maximum binding, about 15 MeV, is 
realized at r s ~ 0.4. One should bear in mind, that in the case of baryon-rich matter local 
minima of e result from a strong cancellation between the attractive scalar and repulsive 
vector interactions. Therefore, it is very sensitive to their relative strengths. The results 
presented above are obtained for Gy = 0.5 Gs- For comparison in Figs. 4 and 5 we also 



The analogous conclusion has been made in Ref. [17] 
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present the model predictions for Gy = 0. In this case the maximum binding energy 
increases to about 90 MeV per baryon and the corresponding r s value shifts to about 0.6. 
It is interesting to note that for Gy = the bound state appears even in the pure u, d 
matter. The corresponding binding energy is about 20 MeV per baryon. 

The dots in Fig. 5 indicate the positions of some conventional baryons. By inspecting 
the figure, one can make a few interesting observations. First, conventional baryons are 
more bound than strangelets even at Gy = . This is an indication that baryon-like 
3-quark correlations might be indeed very important in the baryon-rich quark matter. 
Second, the bound state energies grow monotonously with r s . 



IV. QUARK MATTER AT FINITE TEMPERATURES 

In this section we study properties of deconfined matter at finite temperatures. The 
calculations for this case can be done by using general formulas of Sect. II with the quark 
and antiquark occupation numbers given by Eqs. flOP-fllPP. Unless stated otherwise, the 
results given below correspond to the ratio Gy/Gs = 0.5 . 

Figs. 6(a) represents the pressure isotherms for the case of zero net strangeness 
(r s = 0), which is appropriate for fast processes where net strangeness is conserved, e.g. in 
relativistic nuclear collisions. One can see a region of spinodal instability, d p P < 0, which 
is characteristic for a first order phase transition. The corresponding critical temperature 
is about 35 MeV. The dashed line (binodal) shows the boundary of mixed phase states. 
Bound (zero-pressure) states exist only at temperatures below 15 MeV. 

Generally, the equation of state of the chemically equilibrated quark matter is char- 
acterized by two quantities: net baryon charge, B, and net strangeness, 5*. Therefore, it 
is interesting to study thermal properties of this matter at S ^ . Such states can be 
reached in neutron stars. They can also be realized via the distillation mechanism accom- 
panying a QCD phase transition in heavy-ion collisions p8| |. Fig. 6(b) shows the pressure 



isotherms for r s = 0.4. As discussed above, this case corresponds to most bound strange 
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matter at T = . At this value of r s bound states exist at T < 30 MeV. The dashed 
line in Fig. 6(b) again shows the boundary of the two-phase region. It is obtained by the 
solution of the Gibbs conditions = P^ 2 \ //W — ^(2) ; ^00 — ^( 2 ) where indices 1 and 2 
denote two coexisting phases. As compared to the case r s = the chiral phase transition 
takes place in a wider region of T and ps ■ Applying the Gibbs conditions, one can see 
that for r s = the equilibrium pressure isotherms are constant in the mixed phase region 
(the Maxwell construction). But this is not the case for r s = 0.4. In accordance with 
general conclusions of Ref. p5 |, in the case of two conserving charges (baryon number 
and strangeness) the Maxwell construction is modified in such a way that the equilibrium 
pressure at fixed T increases with pB in the mixed-phase domain. However, this increase 
is small (several percents) and hardly visible in Fig. 6(b). For example, for T = 30 MeV 
the equilibrium pressure changes from 5.36 to 5.44 MeV/fm 3 . It is interesting, that local 
values of r s in the coexisting phases are slightly different (r^ < 0.4 < r^) and only the 
global strangeness ratio is fixed (S/3B = 0.4). 

Fig. 7 shows the critical temperatures for the existence of phase transitions and bound 
states in the equilibrated strange matter as functions of r s . One can see that both 
temperatures first grow with r s and then drop to zero at r s ~ 0.8. The maximal values 
of respectively 50 MeV and 30 MeV are realized at some intermediate r s ~ 0.4 . As 
demonstrated earlier in Fig. 4, this value of r s corresponds to the most bound state of 
strange quark matter at T = 0. So we see an obvious correlation: the deeper is a bound 
state at T = the stronger is a phase transition at finite temperatures. 

It should be emphasized here again that the thermal properties of asymmetric baryon- 
rich quark matter are very sensitive to the relative strength of scalar and vector interac- 
tions |i~9| . If we take Gy = 0, as in most calculations in the literature, the corresponding 
critical temperature at r s = increases to about 70 MeV. On the other hand, if one takes 
Gy = 0.65 Gs, the zero-pressure states disappear completely. The calculation shows that 
there is no phase transition at Gy > 0.71 Gs- It is interesting to note that in all cases 
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this first order phase transition occupies the region of densities around the normal nuclear 
density p • 

A more detailed information about the first order (chiral) phase transition predicted 
by our model is given in Figs. 8 where the critical temperature T c and baryon chemical 
potential \xb = p« + 2//^ = 3/x u are shown for different values of r s . Again one can see 
that maximal T c ~ 50MeV corresponds to r s ~ 0.4 . 

Two phases coexisting in the mixed-phase domain of the chiral phase transition are 
characterized by different values of constituent quark masses. If qq matter evolves from a 
high to a low density state, crossing the two-phase domain, the transition from chirally 
restored states (with low quark masses mf ~ n^of) to the chirally broken phase (where 
rrif ~ m™ c ) takes place. Therefore, a part of the internal energy should be transformed 
into the rest mass. If the total entropy Q S and the effective number of degrees of freedom 
are conserved, the temperature will decrease during this transition. This is clearly seen 
in Fig. 9, where isentropes with entropy per baryon S/B = 2 and S/B = 5 are shown 
for the case r s = . The analogous effect has been predicted within the linear a model 
in Ref. |30| . The drop of temperature as a result of this cooling mechanism may serve 
as an observable signature of the chiral phase transition in nuclear collisions. The above 
behaviour is qualitatively different from the one expected in the deconfinement transition 
where the number of degrees of freedom is smaller in the dilute (hadronic) phase and the 
temperature may increase during hadronization ||31|| . 

These results demonstrate that the chiral phase transition is rather similar to a liquid- 
gas phase transition in normal nuclear matter. The critical temperature and baryon 
density in the present case (T c ~ 30 MeV, pb c ~ Po) are n °t so far from the values 
predicted by the conventional nuclear models j32| (T c ~ 20 MeV, p# c ~ 0.5 p ) . One may 
expect that in a more realistic approach, taking into account nucleonic correlations, the 



5 This quantity should not be confused with the net srangeness introduced in Sec. II. 
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chiral transition may turn into the ordinary "liquid-gas" phase transition. If this would be 
the case, one should be doubtful about the possibility of any other QCD phase transition 
of the liquid-gas type at a higher baryon density. At least only one phase transition of 
this type is predicted within the NJL model. 

V. DISCUSSION OF DECAY MODES 

Let us discuss briefly possible decay channels of bound states of quark matter de- 
scribed above. In strange quark matter (without antiquarks), flavour conversion is only 
possible through weak decays. As follows from Fig. 3, at densities corresponding to zero 
pressure, the condition /i s > /i u holds. This means that weak processes of the types 
s — > u + e~ +V e , s + u — > u + d are allowed. Since there is no local barrier at any r s (see 
Fig. 5), in a system produced initially at some r s ^ 0, all strange quarks will eventu- 
ally be converted into light u, d quarks. Schematicly this conversion process is shown in 
Figs. 9(a) where the initial (r s ^ 0) and final (r s = 0) states corresponds to the left and 
right diagram, respectively. 

This picture is very different as compared to the one based on the MIT bag model. 
Because the quark masses are kept constant (mj = mo/), the condition fi s > m s will 
be first satisfied at a relatively low baryon density oc m^ s . At higher densities a certain 
fraction of s quarks will always be present in a /3-equilibrated matter (see Fig. 9(b)). 

In contrast, in the NJL model the s quark mass is a function of both baryon density 
and strangeness content. As one can see from Fig. 9(a), at any given r s the condition 
fi u ,d = Us can be satisfied only at sufficiently high baryon densities which correspond 
to positive pressure. On the other hand, at the points of zero pressure we always have 
fi u ,d < A* s - Therefore, weak decays will proceed until a system reaches r s = (see right 
picture in Fig. 9(a)). 

The life times of strangelets at T = can be roughly estimated by analogy to the 
neutron decay. The matrix element of s-quark (3 decay is proportional to (A_E) 5 / 2 sin 6 C 
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where AE is the energy gain in the reaction and 6 C is the Cabibbo angle (sin 6 C ~ 0.22 |33| ) . 
As follows from our calculations (see Fig. 3), the energy gain AE = fi s — \x u depends 
on r s , ranging from about 200 MeV at r s = 0.4 to 100 MeV at r s = . By scaling 
with corresponding quantities for the neutron decay, one can express the life time of a 
strangelet as 

/ Am \ 5 

T ~ Tn \AE) Sin ~^ c ' (39) 



where Am = m n — m p ~ 1.2 MeV and r n ~ 880 s is the life time of neutron [33]. This 
estimate gives r ~ 10~ 9 -r- 10~ 7 depending on r s . It is not surprising that this is close to 
the life times of charged pions and conventional hyperons. 

VI. CONCLUSIONS 

By using the NJL model, we have investigated the equation of state of chemically 
equilibrated deconfmed quark matter at various temperatures, baryon densities and 
strangeness contents. The model predicts the existence of loosely bound, negatively- 
charged strangelets with maximal binding energies of about 15 MeV per baryon at 



r s ~ 0.4. Similarly to Ref. [JL7]], no absolutely stable strange quark matter has been 
found. The estimated life-times of these states may be as long as 10~ 7 s. It is shown that 
properties of baryon-rich quark matter are very sensitive to the relative magnitude of 
the vector and scalar interactions. At the standard values of vector and scalar couplings, 
Gv/Gs = 0.5, the metastable bound states of chemically equilibrated matter exist at 
T < 15 MeV, while at Gy = this temperature increases to 40 MeV. 

The calculations reveal the first order chiral phase transition at finite baryon densities 
and moderate temperatures. In the case of zero net strangeness (r s = 0) the critical 
temperature is in the range of 30 MeV and the critical baryon density is around po- We 
believe that this phase transition is reminiscent of the ordinary liquid-gas phase transition 
in nuclear matter. We have found that the maximal critical temperature T c ~ 50 MeV is 
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reached for the ratio of net strangeness to baryon charge S/B ~ 1.2 . The model predicts 
a strong cooling of matter during the phase transition due to generation of the constituent 
mass. 

We hope that these results will generate a certain optimism in searching for unusual 
states of matter in relativistic heavy-ion collisions. 
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1 2 3 4 5 6 



Pb/Po 

FIG. 1. Energy per particle e in pure quark matter at zero temperature as function of 
baryon density at different values of strangeness fraction r s . Points indicate local minima of 
e. po = 0.17 fm~ 3 is normal nuclear density. 
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FIG. 4. Binding energies per baryon in pure quark matter as functions of strangeness frac- 
tion r s . Dotted line shows the results of calculations when the vector interaction is switched 
off (G v = 0). 
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FIG. 5. Minimal energies per baryon in pure quark matter as functions of strangeness frac- 
tion r s . Dashed line shows the same energy in the limit of zero particle densities, Eq. (|36|). 
Different parts of the solid line correspond to metastable (AB and CD) or bound (BC) states. 
Filled squares show masses of lightest baryons. Open squares represent the spin-isospin averaged 
masses of these baryons, e.g. < N >= [mjv + 4 m, a( 1232 )]/5 (for details see Ref. [19]). Dotted 
line in the lower plot shows the results in the limit Gy — ► 0. 
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FIG. 6. Pressure isotherms for chemically equilibrated quark matter with r s = (upper part) 
and r s = 0.4 (lower part). Temperatures are given in MeV near the corresponding curves. 
Boundaries of spinodal regions are shown by the dashed lines. The dashed-dotted line in lower 
plot shows the equilibrium pressure in the mixed phase-domain at T = 30 MeV. 
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FIG. 7. Critical temperatures for existence of bound states (P < 0) and phase transition 
(d p P < 0) in equilibrium quark matter as functions of strangeness fraction. 
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FIG. 8. Critical temperatures and baryon chemical potentials in equilibrated quark matter at 
fixed strangeness fractions r s (indicated by figures near corresponding points). 
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Pb/Po 



FIG. 9. Isentropes S/B = 2 (solid line) and S/B = 5 (dashed-dotted line) of quark matter 
with zero net strangeness. Lower (upper) part corresponds to the /ie — T (p^ — T) plane. Dotted 
and dashed lines represent, respectively, spinodal and binodal boundaries of the two-phase region 
(shown by shading in the lower plot). Point C marks the critical point. 
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(a) 




FIG. 10. Schematic pictures of energy levels (shown by shading) occupied by light and strange 
quarks in cold quark matter at different strangeness fractions r s . Left and right upper di- 
agrams (a) shows the results for strange (r s = 0.4) and nonstrange (r s = 0) matter predicted 
within the NJL model. Low panel (b) represents the same results within the MIT bag model [19]. 
Hatched boxes in upper and low panels shows, respectively, the constituent and bare quark 
masses. Arrows show weak decay processes s — ► u + e + V e and s + u — > u + d. 
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